Heisenberg antiferromagnet on the square lattice for S > 1 
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Theoretical predictions of a semiclassical method - the pure-quantum self-consistent harmonic 
approximation - for the correlation length and staggered susceptibility of the Heisenberg antiferro- 
magnet on the square lattice (2DQHAF) agree very well with recent quantum Monte Carlo data 
for S = 1, as well as with experimental data for the S = 5/2 compounds Rb2MnF4 and KFeF4. 
The theory is parameter-free and can be used to estimate the exchange coupling: for KFeF4 we find 
J = 2.33 ± 0.33 meV, matching with previous determinations. On this basis, the adequacy of the 
quantum nonlinear a model approach in describing the 2DQHAF when S > 1 is discussed. 



We consider the Hamiltonian 
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where i = (ii,? 2 ) runs over the sites of a square lattice, 
i + d defines the four nearest neighbours (ii ± 1, i% ± 1) 
and the quantum operators S\ obey the angular momen- 
tum commutation relations [S'PjSj 3 ] = iS^S^e a ^' r with 

\Si\ 2 = S(S + 1). When J > Eq. @) describes the 2D 
quantum Heisenberg antiferromagnet (QHAF). 

From the experimental point of view the interest in 
this model is due to the existence of several compounds 
that, despite their complex structure, can be described 
by the Hamiltonian (Q) as far as their magnetic behaviour 
is concernedu. Some of them are parent compounds 
of high-T c superconductors, which makes the study of 
their thermodynamic properties of particular relevance. 
Another interesting feature of this class of materials is 
that it contains compounds with 5=1/2 (La2Cu04, 
Sr 2 Cu0 2 Cl 2 ), S = 1 (La 2 Ni0 4 , K 2 NiF 4 ) and S = 5/2 
(KFeF4, Rb 2 MnF 4 ) thus allowing an experimental anal- 
ysis of the 2DQHAF as its quanticity varies from the 
extreme quantum case S = 1/2 to the almost classi- 
cal one S = 5/2. This is of particular interest when a 
comparison with theoretical results is attempted, as the 
spin value S appears as a parameter that can be easily 
changed in most theoretical approaches, obviously save 
the numerical simulations. 

The theory of the QHAF has been generally related 
with that of the quantum nonlinear a model (QNLerM) 
by Chakravarty, Halperin and Nelson (CHN)cl whose 
work led to the first direct comparison between exper- 
imental data and the results of the QNLerM field theory; 
the surprisingly good agreement found caused an intense 
activity, both theoretical and experimental, about the 
subject. In the last few years, however, it turned out 
that for larger spin neither the CHN formulas nor the 



improved ones derived by Hasenfratz and Niedermeier 
(HN)El, can reproduce the experimental data. The dis- 
crepancies observed may be due to the fact that the real 
compounds do not behave like 2DQHAF or to an ac- 
tual inadequacy of the theory. In particular the CHN- 
HN scheme introduces two possible reasons for such in- 
adequacy to occur: the physics of the 2D QHAF is not 
properly described by that of the 2DQNL<tM and/or 
the two(three)-loop renormalization-group expressions 
derived by CHN(HN) do hold at temperatures lower than 
those experimentally accessible. 

The situation can be clarified by using an indepen- 
dent theoretical method, directly applicable to the QHAF 
and whose validity in the temperature region of inter- 
est can be checked. The high-temperature expansion 
(HTE) technique is characterized by such requisites and 
the first well sound doubts about the QNLerM picture of 
the QHAF indeed arose from HTE resultsQ; however, the 
HTE is not applicable in the whole temperature range 
where data from experiments or QMC simulations are 
available, so that it cannot be used to develop a complete 
analysis of the subject. On the other hand, such analysis 
can be carried out by means of the pure-quantum self- 
consistent harmonic approximation (PQSCHA)tJ, whose 
results for S > 1 are fully reliable at all temperatures 
(except the extremely low ones in the 5 = 1 case), as we 
show below. 

The PQSCHA is based on the path-integral formula- 
tion of quantum statistical mechanics, and has been suc- 
cessfully applied to many magnetic systemsO; it permits 
to express quantum thermal averages of physical observ- 
ables in the classical-like form of phase-space integrals, 
where the integrand functions, depending on both T and 
S, are determined from the quantum operators according 
to a precise procedure. The final formulas for the ther- 
modynamic quantities do not contain parameters other 
than those appearing in the Hamiltonian of the model 
under investigation. In the case of magnetic systems the 
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method applies directly to the spin model, with the spin 
value 5 appearing as a coupling parameter that can be 
easily varied. The PQSCHA expression for the statis- 
tical average of a physical observable described by the 
quantum operator O is given by the phase-space integral 
(6) = (1/Z) J d N s O exp(-/3?4ff) where (3 — T -1 ; N 
is the number of lattice sites, s is a classical vector on 
the unitary sphere [\s\ = 1) and Z= J d N s exp(— f3H. c s ) 
is the partition function. The determination of the effec- 
tive Hamiltonian Tl c f{ = 7i c ff({*i}) represents the core 
of the application of the PQSCHA method. The func- 
tion O — 0({si}) is obtained starting from the quantum 
operator O and following the same procedure used to 
determine the configurational part of H. c ff . 

The energy scale JS 2 with 5 = 5 + 1/2 naturally 
appears in deriving the effective Hamiltonian, and we 
hence define, and hereafter use, the reduced tempera- 
ture t = T/JS 2 . In the specific case of the 2DQHAF 
described by Eq. ([!]) we find 
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where 2 = 1 - X>/2, V = £ k (l - 7 £) 1 /2£ k /(jV5) ) C k = 
coth/ k - fe 1 , / k = Wk/(25t), 7k = (cosfci + cosfc 2 )/2 
and k = (fci,^) wave vector in the first Brillouin zone; 
Q (t) is a uniform term that does not affect the evaluation 
of statistical averages. The self-consistent solution of the 
two coupled equations Wk = 4k 2 (1 — 7 2 .) 1 / 2 and n 2 = 6 2 — 
t/(2K 2 ) gives us all the ingredients needed to evaluate the 
thermodynamic properties of the system. Details about 
the derivation of the effective Hamiltonian and of the 
above formulas are given in Ref. § 

The renormalization coefficient T> = X>(5, t) measures 
the strength of the pure-quantum fluctuations, whose 
contribution to the thermodynamics of the system is the 
only approximated one in the PQSCHA scheme: The the- 
ory is hence quantitatively meaningful as far as T> is small 
enough to justify the self-consistent harmonic treatment 
of the pure-quantum effects. In particular the simple cri- 
terion V < 0.5 is a reasonable one to check the validity 
of the final results. 

In Fig. |l| we show the coefficient X>(5, t) as a function 
of temperature and for different spin values: Besides the 
obvious observation that the temperature range where 
V < 0.5 depends on the spin value, we also note that for 
5 = 1 such interval extends to almost the whole temper- 
ature range, leaving the extreme quantum case 5=1/2 
the only delicate one as far as the validity of the PQSCHA 
is concerned. Fig. |l| should clarify that the PQSCHA is 
not a high-temperature method, but rather a scmiclassi- 
cal one whose results are fully reliable already for 5 = 1 . 

The 5 = 1/2 case is extensively discussed in Ref. |?] 
and we do not deal with it in this paper; however, we 
note that our 5=1/2 results agree with QMC data for 
t > 0.4 which means, as seen in Fig. |, V < 0.47. This 
confirms the criterion adopted to be well sound. 




FIG. 1. Renormalization coefficient T> vs t for S — 1/2, 1, 
5/2 (from the top curve). 

From Eq. (||) we see that the symmetry of the Hamil- 
tonian is left unchanged so that the quantum system 
essentially behaves, at an actual temperature t, as its 
classical counterpart does at an effective temperature 
^cff = t/Q^if) . Being 9 4 (t) < 1, we see that the sys- 
tem is more disordered than its classical counterpart, as 
the pure-quantum fluctuations make the former behave 
like the latter does at a higher temperature. 

As for the correlation functions G(r) = (SfS %+,.), with 
r = (r%,r2) any vector on the square lattice, we find 
G(r) = 5 2 ^(si-Si +r ) c ff where ( • ) e fi is the classical-like 
statistical average with the effective Hamiltonian; the pa- 
rameter 9 4 = 1 — XV, defined in Ref. ||, goes to a constant 
and finite value for large |r|. From the above formulas, 
the correlation length, defined from the asymptotic ex- 
pression G(r) oc exp(— |r|/£) for large |r|, turns out to 
be 



m = &l(teff) 
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where £ c i is the correlation length of the classical HAF, 
unique ingredient we need to obtain £ (t) for the quantum 
system, being the evaluation of 9^(t) a simple matter for 
any spin value. The PQSCHA expression for the stag- 
gered susceptibility \ = S r (~ ) ri+r ' 2 G(r)/3 is 
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and in this case we need to know the classical G(r) at 
any r to obtain the numerical value of the quantum \. 

Figs. 11^1 show our results together with the available 
QMC and experimental data. We underline that no 
best-fit procedure is involved in such comparison as the 
PQSCHA has no free parameters once J and 5 are given. 

In the 5 = 1 case (Fig. ||) we compare our curves for 
£ and, x with the new QMC data obtained by Harada 
et al.a; such data, which unfortunately do no extend to 
low temperatures, do in fact sit on our curves. Also the 
experimental data for La2Ni04 and K2N1F4, which are 
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not included in Fig. g for the sake of clarity, very well 
agree with our PQSCHA curves as shown in Ref. [|. 
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FIG. 2. Correlation length £ and staggered susceptibility 
X* = x/S 2 vs t for S = 1. Symbols are QMC data from 
Ref. |. 

For S 5/2 we use the experimental data relative 
to KFeF4Ll and also those for Rb2MnF4 very recently 
obtained by Lee et al.Eii Both compounds are char- 
acterized by an Ising-like anisotropy. H A ~ J/20: Fol- 
lowing the reasoning by BirgeneauEll, the crossover be- 
tween Ising and Heisenberg behaviour should occur when 
h A £ 2 ~ 1, with h A = H A /{zJS) ~ 0.005. These ma- 
terials are hence expected to behave like 2DQHAF for 
£ < 14, i.e. t > 0.6, which is in fact the region where 
the PQSCHA curves agree with the experimental data, 
as seen in Figs. [|^. 

As for the compound KFeF4 a few more comments are 
in order; its magnetic ions are distributed on a non per- 
fect square lattice and the magnetic interaction is hence 
characterized by two different exchange integrals J a and 
Jb; such difference is small enough to allow the system 
to be described by Eq. (|l|) but with the value of J, as 
from different experiments, slightly variable. In particu- 
lar from susceptibility, Raman scattering, and neutron 
scattering measurements J is found to be 2.30 meV, 
2.35 meV, and 2.44 meV, respectively. Neutron scat- 
tering experiments by Fulton et alrl measured J a and Jb 
separately getting J a = 2.18 meV and Jb = 2.73 meV at 
50 K and J a = J b = 2.40 meV at 100 K. In our theory, 
on the other hand, there is no free parameter and the 
exchange integral J only enters the definition of the re- 
duced temperature; its value can be hence easily derived 
by optimizing the agreement between the experimental 
data for the correlation length and our curve. In this 
procedure, in order to take into account the above men- 
tioned effects of the Ising anisotropy, we have only used 
the experimental data with £ < 14, obtaining the value 
J = 2.33 ±0.03 meV; this result does not change discard- 
ing the data with 8 < £ < 14 and agrees with the above 
mentioned independent determinations. 




FIG. 3. Correlation length £ vs t for S ff 5/2. Symbols 
are experimental data for KFeF4 (circles)u and Rb2MnF4 
(squares)E3. The inset shows a magnification of the region 
where the crossover between Heisenberg and Ising behaviour 
is observed. 
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FIG. 4. Staggered susceptibility = x/S 2 vs t for 
S — 5/2. Symbols and inset as in Fig. 

From the extensive comparison between our curves and 
all the available experimental and QMC data, as well as 
with the HTE results as shown in Ref. ^, it is evident 
that the PQSCHA leads to a proper description of the 
thermodynamic behaviour of the QHAF. 

Let us now consider the temperature dependence of the 
correlation length in the S = 1 case (a similar analysis 
can be developed for any S > 1): Our data extend up 



to £ 
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i.e. t 



0.2, where the renormalization 



coefficient T>(S = 1) is still less than 0.45. This means 
that we can safely assume our curves to reproduce the 
correct £ for any t > 0.2. Notice that this lower limit is 
set by the absence of classical values for £ > 10 5 . 

Having said that, we try to fit our curves with the 
low-temperature CHN-HN formula 
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where p and c are the two fitting parameters that do not 
depend upon T. In order to assert that Eq. (|j) describes 
the correct temperature dependence of £ in the tempera- 
ture range of interest, the fit must be stable, in the sense 
that the resulting values of p and c must not vary if the 
fit is restricted to lower or higher temperatures. In what 
follows we show that, in fact, this is not the case. 

We have fitted our curve with Eq. (^) in different in- 
tervals of temperature and found that the values of p and 
c change drastically when the low temperature data are 
not included in the fit. To enlighten the discrepancies in 
Fig. H we report the quantity t ln£ as a function of t; the 
full line is the PQSCHA curve while the other curves are 
given by Eq. (|]) with p and c as from the fit in the five 
temperature intervals [x, 1.7] with x = 0.2, 0.46 and 0.56. 
The inset shows p* = 2np/S 2 , as resulting from fits in 
the temperature intervals [x, 1-7], as a function of x. The 
fit is clearly unstable and Eq. (Q) cannot reproduce the 
correct behaviour of the 2DQHAF correlation length in 
the whole temperature range. 
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FIG. 5. t\n£ vs t for S = 1: fits of the PQSCHA curve 
(full line) with Eq. (^]) in the temperature interval [x, 1.7] with 
x — 0.2 (dashed line), 0.46 (dotted line) and 0.56 (dash-dotted 
line). The inset shows p = 2ivp/S 2 as from fits in the temper- 
ature interval [x, 1.7] as a function of x. 

From this analysis we conclude that the CHN-HN the- 
ory does not describe the QHAF when £ < 10 5 , i.e. for 
t > 0.2 for S = 1 and, more in general, for t > to where 
to is such that £ c i(W^(*o)) ~ 10 5 for any S > 1; the only 
reason why Eq. (0) seems to reproduce the experimental 
data is that the temperature range where the fit is carried 
out is small enough not to show significant deviations. 

The CHN-HN theory is in fact already known to de- 
scribe the QHAF only at sufficiently low temperatures, 
but no quantitative indication is given about the actual 
range of validity, save the fact that this is narrower for 
larger spin. What our work clearly shows is that such 
range of validity, for 5 > 1 , does not overlap with the re- 
gion where experimental and QMC data are available; in 
other terms, such data are not described by the three- 
loop renormalization-group expressions relative to the 



QNLctM. 

Although our analysis cannot be extended to all tem- 
peratures in the 5 = 1/2 case, we point out that con- 
clusions similar to ours have been recently drawn by 
Beard et al.E£l from their low-temperature QMC results. 
On the other hand, Kim and Troyer, who also per- 
formed QMC simulations on the 2D QHAF with 5=1/2, 
clain£3 that their results are in "excellent agreement" 
with the QNLoM field-theory predictions. Nevertheless, 
such agreement involves unstable fitting procedures; in 
particular, the uniform susceptibility Xu can be fitted to 
the CHN-HN expression only for t < 0.23 while for the 
correlation length the restriction £ > 39.2 (i.e. t < 0.27) 
is necessary to make the fit stable but not even sufficient 
to let the resulting values of p and c coincide with those 
obtained by fitting \ u . In fact we think that the results 
presented by Kim and Troyer, which moreover extend to 
temperatures not as low as those studied by Beard et 
al., do not "confirm the validity of the mapping from the 
QHAF to the NLcrM", but rather suggest that also for 
5 = 1/2 the CHN-HN formulas do not properly describe 
the behaviour of the 2D QHAF for t > 0.25, i.e., in the 
temperature region where most experimental and QMC 
data are available. 

We are grateful to Prof. R.J. Birgeneau and to Y.S. 
Lee (MIT) for fruitful exchanges and for providing us 
with experimental data prior to publication. 
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